unvariance <- function(data, na.rm = TRUE) {
	return(sd(data,na.rm=na.rm)^2)
}

